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Abstract 

We develop a semiclassical method to calculate the density of magnetic 
monopoles in non-abelian gauge theories at finite temperature in the dilute gas 
approximation. This quantity is related to the inverse magnetic screening length 
for which we obtain fXM = 0.255g^T in SU(2). 



1. Introduction 

The screening of static magnetic gauge fields in Yang-Mills theories at high temper- 
ature has long proved to be intractable by analytic techniques. There exist well known 
arguments^ that the inverse screening length hm must be of order g^T, and this has been 
partially supported by numerical calculations in the framework of lattice gauge theory^ 
which have yielded the result fXM = (0.27 ± 0.03)g^T for SU(2). The recently developed 
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resummation techniques for the finite-temperature gauge theories cure many infrared di- 
vergences arising in the perturbative expansion, but not those associated with long-range 
static magnetic gauge fields. Although the Schwinger-Dyson equation that determines /im 
has been identified some time ago^, its solution remains unknown. A recent attempt^ to 
calculate magnetic screening to leading order in dual QCD has led to the result 2nT/g 
for the dual gluon mass in the high-temperature limit, which is not directly related to the 
static magnetic screening mass. 

Here we propose a new approach to the calculation of fiM, which is not based on a 
perturbative expansion in the gauge coupling constant g. Our starting point is the obser- 
vation* that the combination g'^T defines an inverse length without any factor involving 
powers of Planck's constant h. One may therefore speculate that //m can be calculated 
from classical Yang-Mills statistical mechanics, with quantum effects providing corrections 
of higher order in g. This conjecture was recently found to be true for the thermal gauge 
boson (plasmon) damping rate, which also is of order g'^T and can be calculated either by 
resummation techniques^ or by a stability analysis of time-dependent classical gauge field 
configurations^. This coincidence is nontrivial, because in the perturbative calculation the 
gluon damping rate is seen to superficially depend on the electric screening mass, which is 
of order gT/V^- However, this dependence exactly cancels in the leading order result. 

Similarly, it has been emphasized by Landsman® that the effective high-temperature 
effective action for the static sector of gauge theories contains quantum corrections describ- 
ing the Debye screening of static electric fields. These effectively convert the dimensionally 
reduced Yang- Mills action into a Yang-Mills-Higgs action that remains unbroken at tree 
level. Magnetic monopoles can induce a symmetry breaking term into the effective Higgs 
action, as pointed out by Oleszczuk and Polonyi^. Our treatment assumes that these 
* Note that in our convention g is defined by the classical Yang- Mills equations, and is 
related to the dimensionless coupling constant by ag = g'^h/Air. 
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effects of the full quantum theory do not contribute to the magnetic screening mass at 
leading order in g. Whether that is indeed so will require futher investigation. 

This paper is organized as follows: we connect the magnetic screening mass in high- 
temperature Yang-Mills theories with the magnetic monopole charge density in the frame- 
work of the linear response theory. We then calculate the density of magnetic monopole 
charges from the canonical partition function in Chapter 3. A saddle point evaluation in 
the space of general radial symmetric monopole field configurations becomes possible by 
inserting an appropriate unit factor, which sets the scale and stabilizes classical soliton 
solutions of a given size. The integration over scale sizes is done at the end. Besides these 
stabilized classical field configurations the contributions of fiuctuation modes with zero 
frequency, i.e. the collective modes of translation and rotation, must be taken into account 
because they determine the entropy of a given field configuration. Finally, in Chapter 4, 
we conclude by exploring the characteristic size and energy of dominant monopole config- 
urations under conditions that may be reached in heavy-ion collisions at RHIC energies. 

2. Magnetic Screening 

In order to obtain the magnetic screening length one usually investigates the linear re- 
sponse of the medium to an infinitesimal external magnetic field. The magnetic field, which 
couples to the monopole charge, is curl-free and can be expressed by a scalar magnetic 
potential 0", 

Bt = -dir- (1) 
The application of such an external magnetic field modifies the monopole charge density 
in the plasma, so the partition function becomes 

Z(0) = TVe-^(^-'^"^"), (2) 

where is the magnetic charge operator and the trace "Tr" runs over all possible states 
of the Yang- Mills field. Here we discuss how to evaluate this partition function in the case 
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of SU(2) only, but the result can easily be generalized to SU(N). Let us denote a given 
color charge state by \q, q^), where q, the eigenvalue of the Casimir operator, is a multiplet 
index. A magnetic monopole state corresponds to the eigenvalue q = 1- 

Since the Yang-Mills action is color symmetric the eigenstates of the Hamiltonian are 
gs independent, so one can always "rotate" an energy and color eigenstate \u), q^) so that 
the action of (3(j)°-Q°- can be represented by that of the generator of the abelian subgroup 
©/3 = (7~^/3|0|/3, where is the elementary magnetic charge, has the eigenvalues 
= — 1, and +1 and finally 

101' = (3) 

Using this particular representation of gluon states 1^ and H commute and in the dilute 
gas approximation the quantum numbers uj and q^ are additive, so the partition function 
(2) factorizes for different uj and gs values. -"^"^ Taking into account the bosonic nature of 
excitations of the Yang- Mills field, the magnetic monopole partition function in the dilute 
gas limit becomes: 

oo 

where q^ runs over the possibilities —1, and 

The summation over the indefinite occupation number n yields a Bose-Einstein factor 
for each uj and 53 . From this partition function we obtain the magnetic color charge density 
using the standard formula 

1 d 

\nZ. (5) 
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In the weak external field limit (0 <S 1) this leads to (see Appendix A for details): 
for the group SU(N) with the effective one-monopole partition function 

UJ 



Because of the large average mass of a monopole we can make use of the Boltzmann 
approximation, obtaining 



UJ 



Y^e-f^ = Tr{e-^^). (8) 



Since the magnetic monopole charge density is the source of the divergence of the 
magnetic field 

diB'i = p«, (9) 

we arrive at the following equation describing the linear polarizability of a magnetic 
monopole gas 

did^r - -^r = 0. (10) 

Inspecting this equation one easily realizes that an effective magnetic screening mass is 
obtained, 

m=(T^)^9'T, (11) 



causing an exponential damping of the external magnetic potential 

0«(r) ~ e-'^^^ (12) 

It is reasonable to expect that the ratio Zm/V scales with at high temperatures 
because T is the only available scale. Recall furthermore that the four-dimensional gauge 
theory at high temperatures undergoes dimensional reduction and becomes equivalent to 
a three-dimensional gauge theory with the effective temperature g^T. We therefore expect 
that the partition function scales as Zm ~ Vig'^T)^, so the magnetic screening mass is 
proportional to g'^T. 

It is illuminating to repeat the same considerations for the electric sector of the gauge 
theory. In this case the external potential must be taken as that of a static color-electric 
field, and one considers the dilute gas of gauge field excitations with electric charge in the 
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adjoint representation. In lowest perturbative order these are just the free perturbative 
gauge bosons, but it is now inappropriate to use the Boltzmann approximation because 
there is no mass gap. As shown in Appendix A, this consideration yields precisely the 



3. Partition Function of Magnetic Monopoles 

3.1. Spherically Symmetric Gauge Fields 

In order to calculate the nontrivial factor we need to obtain the partition function Zm- 
Although simple scaling arguments^'^^ show that there may be a characteristic monopole 
size Rq contributing dominantly to the partition sum 



no classical monopole solution exists, which would offer a stable stationary point of the 
path integral defining the partition function. 

While the t'Hooft-Polyakov solution^^ for spontaneously broken SU(2) gauge theories 
is stable and has finite energy, its counterpart in pure SU(2) gauge theory, the Wu-Yang 
monopole^"^, has infinite energy and is unstable against small perturbations^^'^^. This lack 
of a basis for a semiclassical analysis of the statistical mechanics of monopole solutions in 
non-abelian gauge theories apparently vitiates an analytical approach. Here we want to 
sketch how the stability problem may be circumvented, opening the way to a semiclassical 
calculation of the magnetic monopole density at high temperature, and therefore of the 
magnetic screening length. 

Let us begin by considering SU(2) gauge field configurations which carry one unit of 
magnetic charge. They must asymptotically look like a monopole {j = 0) mode of the 
operator 



electric screening mass ^^/N/3gT obtained in diagrammatic perturbation theory. 




(13) 



J = S + I + L, 



(14) 
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where S and / denote the generators of spin and color spin in the adjoint representation 
and L is the generator of orbital angular momentum for the gauge field. Since the gauge 
field belongs to the representation S = I = 1, there are three different possibilities to 
construct aj = mode; namely \L,T) = |0,0), and |2,2) combinations, where the 

grand-spin quantum number T is obtained from the eigenvalue of the Casimir operator 

f = S + i. (15) 

The Wu-Yang monopole, known to be unstable, belongs to the j = mode of type \L,T) = 
|1, 1). A further unstable mode has been found to involve a combination of the two other 
states by Akiba, Kikuchi and Yamagida^^. 

Having this in mind we start our investigation with the most general ansatz^^ for the 
monopole vector potential 

A^a ^ E ^a"^^^"^ = l ^(^^''^ ' + (^^"^"^ + ^) + ' (l^) 

a=±,0 

where u{r), 4){r) and w{r) are real functions of the radial variable only and we use the 
projectors 

Pia = riiUa (17) 

with the unit radial vector ni = Xi/r. Configurations of the Wu-Yang monopole-type are 
proportional to Yi^^ia^ -^ia) ~ ^iaj^j: i-G. they belong to the choice w = Q and tt = 0. The 
unstable mode found by Akiba et al. is proportional to a linear combination of |(-Pia + -Pja) 
and Pf^, hence corresponds to w ^ and sin^ = 0. 

The different components of this ansatz, A'^,A~ and A^, the respective coefficients 
of the projectors, are related to the \L,T) specification of the monopole mode {j = 0) by 
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simple linear combinations 



0,0) = 



(18) 



|1,1)=2(^^-^") 



(19) 



and 



|2,2) = i(A++yl-)-i^°. 



(20) 



The magnetic field described by our ansatz is 

1 



B. 



la — ^2 



[P+ {iru' + u{r(j)' + w)) e-"^ + c.c. + (1 - u^)\ 



(21) 



where the prime denotes radial differentiation and c.c. stands for complex conjugate. 

The magnetic monopole charge seen from outside a sphere of radius r can be obtained 
from the magnetic analogue of Gauss' law 



Qa{r) = ^ J d^rdiBia = ^ j^UiBia = na{l - w^(r)). 



(22) 



This result shows that a monopole field configuration requires asymptotically tt — > as 
r ^ oo, irrespective to the fields w{r) and 0(r). 

The energy of a static configuration in the high-T limit defines the effective action of 
the dimensionally reduced euclidean field theory 

1 



(23) 



(3E^(3 dr 
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Introducing the scaled parameter P — we find for the ansatz (16): 



1 



du\ 

-r) + :r + I + 

ar J \ dr r 



2r2 



(24) 



Here one observes that only the field u{r) is really a dynamical degree of freedom 
while w{r) is non-dynamical and 4>{r) is a cyclic variable. Their physical interpretation 
becomes clear inspecting infinitesimal gauge transformations of the vector potential 

6Aia = DiabSAt, (25) 

where 

Diab = Sabdi — eacbAic (26) 

is the gauge-covariant derivative. Any gauge transformation which conserves the monopole 
form of our ansatz (16) must have the form 

6Ab = Ub ■ SA{r). (27) 

This implies that restricting ourselves to static magnetic monopole gauge field configura- 
tions there still remains a residual gauge degree of freedom. The variation of the ansatz 
fields under such an infinitesimal gauge transformation is given by 



Su = 


0, 


5(1) = 


-6A, 


Sw = 


r—SA 




dr 



The field u{r) and hence the monopole charge is unchanged, but the field (p is rotated by 
a gauge transformation mixing fiuctuations of the pure Wu-Yang ansatz with the other 
unstable mode. The meaning of the field w{r) is less obvious, we note only that it vanishes 
in the Schwinger gauge XiAia = 0. 

3.2. The Functional Integral 

Inspecting the form of the effective action S^, one realizes that the path integral in 
the canonical partition sum can be easily done in two of the fields, w{r) and 0(r). Turning 
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from the 'Cartesian' field variables Aia (in mode j = 0) to the 'cylindrical' variables u, (j) 
and w introduced before we use the integration measure 

VA'^^ = ^Vu V(j) Vw. (29) 

Because the ansatz (16) is invariant under transformation by J, the single monopole par- 
tition sum factorizes for small fluctuations around the monopole form: 

i>o 

where 

[u] = j V{A^) exp (-^3 [^^=° + 5A^] + ^3 [^^'^^ ) (31) 

depends only on u{r) due to the gauge freedom (28). Here we have anticipated that we 
will find a nontrivial stationary point in the mode tt(r), around which we can expand the 
functional integration. All integrations in (31) will be Gaussian except those corresponding 
to the zero-modes associated with translations and rotations: 

Z,->o = n [^] = [u] Zrot [u] n 'Zj [u] , (32) 

i>o j>0 

where the prime indicates that the zero-modes have been separated. From studies of small 
fluctuations around the sphaleron^®'^^ in spontaneously broken SU(2) gauge theory it is 
known that all other modes are stable, and lead to well-behaved Gaussian integrals in the 
semiclassical approximation. 

The Gaussian integral in the variable w can be easily done. After doing that we get 

/—Vw exp I —(3 [ dr( — + ^ j | = const.. (33) 
r \ Jo K'T drj J 

compensating for the factor u in (29). We normalize the constant to be unity. Now the 

compact integral over the fleld (f) is trivial and can be normalized such that 

P0=1, (34) 
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resulting in a functional integral over u{r) 



T>u e 



(35) 



if we neglect the influence of non-collective multipole (j > 0) fluctuations on the partition 
function. This issue will be briefly discussed at the end of this Chapter again. 

One may wonder at this point why we have chosen these particular normalizations. 
The physical idea behind it is that for free fields one has to arrive at the canonical partition 
sum at high temperature. Since the only dynamical degree of freedom is represented by 
the other two fields must not contribute to a correctly normalized path integral. 

3.3. Setting the Scale 

Unfortunately, E[u\ does not have a stable saddle point UQ{r). This lack of a basis 
for the semiclassical expansion, however, does not preclude the calculation of the partition 
function for gauge field configurations u with monopole symmetry. It only implies that the 
evaluation of the functional integral (35) cannot be restricted to integration over Gaus- 
sian fluctuations around a classical solution UQ{r). We now propose a method how this 
functional integral (35), where u{r) satisfies the boundary conditions w(oo) = 0, w(0) = 1, 
may be calculated. It may be practically carried out by adding a stabilizing term to the 
expression for E[u]^ 



introduced in eq. (24). The idea is to introduce a length scale that favors monopole 
configurations of a particular core size, and then integrate over the dummy scale in such a 
way that the partition function remains unchanged. Inspired by the analogous expression 
for the t'Hooft-Polyakov monopole^^, where a stabilizing mass term (gf$)^u^ involving the 
Higgs field $ appears, we introduce the term A£'['u] = X^D[u]^ 




(36) 




(37) 
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by inserting a unit factor 

/•OO 

1= / d(A2) /^L'fw] e-'^^'^M (38) 
Jo 

into the integral (35). Here A has the dimension of an inverse length. At large r, in view 
of the boundary condition u{r) '^^^^ 0, it acts like a mass term leading to the asymptotic 
solution u{r) Aexp(— Ar), while it does not interfere with the limit u{0) = ±1. Our 
choice of D[u] is unique* if one requests that the integrand be an even function of not 
higher than fourth order in u which leads to a vanishing energy density at r ^ oo. 

After inserting the unit factor (38) into the functional integral (35) we interchange the 
order of integration over the dummy scale parameter, A, and with that over the fields. A 
consistent use of the dimensionless variable x — \r leads now to the following expression: 



where 



(40) 



dx ) ' x"^ 

D[u] = -j^ dx{2u^-u% (41) 
AE[u] = X'^D[u]. (42) 



* We have explored the consequences of the modified scale breaking factors D[u] = 
J dr{l — (l—u^)^), which are of higher order in u{r). For A/" = 3, 4 these yield somewhat 
smaller values (by 7 and 14 percent, respectively) for the monopole density within our 
approximation. This indicates that our result is not very sensitive to the precise form of the 
scale breaking term. We note that the functional integral over u{r) in (35) could, of course, 
be performed numerically by Monte-Carlo integration, avoiding the errors introduced in 
the saddle-point approximation. 
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The functional integral over u in (39) can now be approximated by a Gaussian integration 
around the lowest energy stationary solution u{r) — uo{x) of the exponent, which satisfies 
the equation 

^^?+Ml-«^J,-l)^0. (43) 

Since there is no other scale involved besides the dummy parameter A, the solution of eq. 
(43) is solely a function of the dimensionless variable x, and the ground state energy scales 
as 

E[uo] + AE[uo] = Xa. (44) 

By numerical integration of (43) we have found the value a = 1.469. The function uo{r) 
is displayed in Figure 1. We also obtain A£^[uo] = bX with b — 0.695. The monopole 
charge contained inside the radius r (eq. 22) is also shown in Fig. 1. The radius R1/2 
inside which half of the asymptotic charge is contained can be defined as the size of the 
monopole-soliton. This happens to be the inflection point of the solution uo{Xr) 

R,/2 = X-\ (45) 

Postponing the Gaussian integration over small fluctuations around uq we are left with a 
single-parameter integral: 

nOO 

Zm-W dXb e-^"^ ZtrZrot. (46) 
^0 

3.4. Zero-Mode Contributions 

Now we turn to the determination of the zero- mode contributions to Zm- The method 
we use is that of collective coordinates, ^° which leads to an integral over the respective 
group volumes of translations and rotations, V and Stt^. What remains is to take into 
account the normalization of the wave functions describing these modes. Physically they 
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are related to the total momentum (energy) of a monopole soliton and to its rotational 
inertia, respectively, as Ztr = X^VM^^., since the classical solution is a function oi x — Ar, 
and Zrot = 87rWrot5 with 

K = ^ld'r 6(r), (47) 

■^rot ~ Q~ J '^^^ — correction) , (48) 

where e(r) is the energy density of the classical soliton, and the correction — due to an extra 
gauge rotation — ensures the correct boundary condition for the rotational energy density 
at radial infinity (see ref. 19). The collective wave function would have the form 

^ = ^ e'^- /, (44) 

with a form factor / normalized so that including a spin factor 3 for massive solitons and 
a radial normalization factor 47r we have 

kl'=-A/'t7'-^-3-47r = l, (50) 

whence we obtain eq.(47). While the total energy of the soliton scales like E = AnXa, the 
rotational inertia of a homogeneous sphere with radius R like ^ER'^ = |47ra/A, because 
the characteristic radius of the soliton is A~^. This yields the approximate expressions 

Af^, ^ ^Xa, (51) 



and 



which lead to the zero-mode factors 



3 

2 ^ 2 



Ztr = A3 Ztr =(30) A^ AV, (53) 



^rot 



3 

2 \ 2 , •3^2 



) A-5 8n\ (54) 
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We note that their product scales hke A^. The integral over A in (46) can now be carried 
out, yielding 

Zm = 12Ztr^rot^7— (55) 
Inserting (53,54) we finally obtain the monopole density in the classical limit as 

PM = = Y^Ag'Tf = 0.0657 ^ (g'Tf. (56) 



Tra a 



3.5. Radial Monopole Fluctuations 

We finally evaluate the scaled determinantal contribution, A~^/^, due to the fiuctua- 
tions around the soliton- monopole solution uo{r): 

A = n^n- (57) 

n 

Here uj^ are the eigenvalues of the operator 

f]^ = -^ + (l-3.g(.))(l-l) (58) 

obtained expanding the effective action E[u\ + AE[u] up to second order in {u — uq). The 
effective potential involved in this eigenvalue problem. 



(59) 



is plotted in Figure 2. Since we are calculating the partition sum Zm, which is restricted 
to single monopole configurations, we only take into account fiuctuations around uq that 
vanish at infinity. Hence only those eigenvalues corresponding to a bound state in this 
effective potential contribute, i.e. 

A = l[u;le{l-u;l) (60) 

n 
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where 9 is the step function. Numerically we have found only one such bound state with 
the eigenvalue uIq — 0.950. 

However, the asymptotic form of V^fp, — is known to have infinitely many bound 
states close to threshold. Their eigenvalues, estimated from the Bohr-Sommerfeld formula. 



are 

ujI=(i- e-(2-+i)-) . (61) 

The spatial extension of the corresponding eigenfunctions (n > 1) is so large, that they 
exceed a radius 70 times larger than the soliton itself. It is therefore questionable whether 
these fluctuations are physically meaningful and should be taken into account at all for 
a monopole gas of finite density. Fortunately, this issue is not really critical, since the 
contribution of such large-size fluctuations is very close to a factor 1 

Y[ujI = exp L ^ 1.0008. (62) 

n=l ^ ^ 

We are left with the contribution of the numerically found ground state in this effective 
potential which — normalizing the same way as before — reads as 

Pi — ^ 1.03. (63) 

Collecting all known factors together we obtain the monopole density 

p = 0.0657 -^{g^Tf « 0.0326 {g'^Tf, (64) 
aujQ 

leading to a screening mass of 

liM ~ 0.255 (g^T) (65) 

for SU(2), which is amazingly close to the value obtained in lattice simulations^. Of course, 
the contributions due to higher multipole modes, which were neglected here (cf. eq. 35), 
may change this result quantitatively. However, the fact that the determinantal factor for 
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radial fluctuations is very close to unity lets us be optimistic that the influence of those 
higher modes on the numerical constant in eq. (65) will be minor. 

4. Conclusions 

Finally we would like to present some numerical estimates characteristic for QCD 
under conditions that may be reached in nuclear collisions at RHIC energies. We assume 
a temperature of T = 300 MeV and a coupling constant of ^ = 2, corresponding to 
ctg = 5'^/47r ~ 0.32. We note that while this value of as would justify a perturbative QCD 
approach, the monopole charge = 0.5 may justify the neglect of monopole-monopole 
interactions. The magnetic screening mass for SU(3) we obtain using the above values is 
liM ~ 375 MeV yielding a screening length of /im = 0.53 fm for static magnetic fields. 

The size and the total energy of an average monopole-soliton can also be estimated 
using the scale representation of eq. (46): 

/•oo 

Pm= / (iAn(A)e-^«^ (66) 
Jo 

with n(A) = uqX^. The average soliton size is obtained as 

RM^{X-')^lap^^, (67) 

which is about 1 fm at this temperature. Unfortunately this exceeds the average dis- 
tance do 4:/g'^T between monopoles at the equilibrium density (64) indicating a possible 
breakdown of the dilute gas approximation employed here. On the other hand, the average 
interaction energy between monopoles 

Eint = (47r^'rfo)-' « 0.02 T (68) 

is quite small compared with the average monopole mass 

{E) = E[uo]/g^ = 4T. (69) 
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This indicates that higher-order terms in the effective action Ss[A] for superpositions of 
monopoles are small, and gives reason to believe that the dilute gas approximation may 
be quite trustworthy. It allows for an intuitive understanding of the mechanism for static 
magnetic screening in thermal gauge theories. We hope to return to the problem of evaluat- 
ing the determinant of eigenvalues of higher modes, as well as to the question of non-static 
loop corrections in the future. 
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Appendix A 

In this appendix we derive the color charge density (eq.6) from the dilute gas partition 
function (4) in the weak external field limit. From eq.(4) we obtain the following expression 
for the logarithm of the partition function 

lnZ = -^^ln(l-e-'5'^e®^3)^ (^j) 

which after factorizing the color charge independent one-gluon partition functions can be 
casted into the form 

lnZ = - J] J] [ln(l-e-'^") + In (l + n(cc;) (l - e®^^))] , (A.2) 

where 

^(^) = (^-3) 
is the Bose-Einstein distribution function. 

Restoring the generality of the discussion we replace now Qqz again by f3(j)°'Q°' and 
write the summation over q^^ as the adjoint color trace tr. We arrive at 

lnZ= - J] [ln(l-e-^'^)trl + tr In (^1 + n(cc;) (^1 - e^'^''^"))] . (A.4) 

Separating now the external field independent term, 

lnZo = -(iV2-l) J]ln(l-e-^'^), (A.5) 

we are left with 

lnZ = lnZo-^trln(^l + n(cj) (l - e^'^''^ ^^^j , (A.Q) 

UJ 

which can be approximated in the dilute limit n{u>) <C 1 by 

lnZf«lnZo-^n(w)tr (^l-e^'^"^") . {A.7) 
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Finally we use the fact that the external field is infinitesimal in case of seeking for a 
linear response of the gluon medium. We obtain 

hiZ ^hiZo + Y^ n{uj) tr (/3(I)''Q'' + y^'^/g^O'' + •••) , (^-8) 

which, normalizing the color magnetic charge generators according to the adjoint repre- 
sentation of the SU(N) algebra, 



leads to 



InZ ^InZo - J]n(a;). {A.IO) 

The definition (5) of the color charge density leads finally to 

p% = —^\nZ ^ ^(/>"Vn(a;), iA.ll) 

which is equivalent to eqs. (6) and (7). 

To calculate the screening of static electric fields, one replaces the scalar magnetic 
potential by a Coulomb potential 

and the magnetic charge operator Q"- by the electric charge operator Q%. The sole change 
in the derivation following equation (A.l) then is that (A. 9) is now replaced by 

tr (q%Q%) = -NgH-K {A.13) 

We finally obtain for the induced color-electric charge density: 
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UJ 



Using the free gluon dispersion relation we have 
hence 

This is the standard perturbative result for the static elctric screening mass in thermal 
Yang-Mills theories. 
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Appendix B 

In this appendix we derive the expressions (21-24) for the magnetic field and energy 
from the ansatz (16) for the vector potential. It is useful here to introduce some short-hand 
notations. We write 

= A+P+ + A-P- + A^Pl. (B.l) 

For calculating derivatives of this expression we start with the projectors. First we note 
that 

djUi = ^{Sij - UiUj) (5.2) 



and 



Sj{nink) = ^{Sij,nk + Sjkrii - 2ninknj). (5.3) 



It is easy to obtain then 



and 



(^ijkdjiekaini) = ^{dia + riiUa) (BA) 



^ijkdjinkHa) = ^eiakUk. (-B.5) 
Using this we obtain the following expressions for the curls of the projectors 

1 Z 1 

"^ijkdjPita = -^^iakTlk + —{^ia + U^Ua) ^ ' {P^a + Pfa) > (^-6) 

^ijkdjPka = -^eiakUk - ^{5ia + rHUa) = ^ {P^~a + Pia) {B.l) 

and 

^rjkdjPl = \^iaknk = {P+ - Pia) ■ (^-S) 

It is now easy to decompose the curl of the vector potential remembering that dj acts on 
pure radial functions as rij . Using 

1 i 

^ijkrijPka = -2^iajnj + -{Sia - riiUa) = iPi^, (5.9) 
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(S.IO) 



and 



we get 



t^jkTljPL = (S.ll) 



= ^P^^ + 1a+ - Ia^^ + C.C + ^Pl - . (5.12) 

We note that = A", (i^^)* = i^;. 

After evaluating the abehan part of the magnetic field (eq. 21) we turn to the calcu- 
lation of the nonabelian part 



(S.13) 



Noting that 



^abc^ijk 



det 



^ai ^aj ^ak 

hi hj hk 

^ci ^cj ^ck 



we obtain 



N,. 



la — 2 



(B.U) 



(S.15) 



Now it is straightforward to evaluate A'^^ and its traces, using the projectors. We get 



and 



{A\, = A+A+P+ + A-A-P- + A^A'P^, 



{A')n^{A+f + {A-f + {Ay. 



(BAG) 
(S.17) 



Inserting these expressions into (B.15) we get 



Nia^6ai{A+A-+A-A^+A^A+)-A+{A-+A'')P+-A-{A++A'')P^,-A''{A++A-)P^,. 

(S.18) 
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Noting that P+ = P-, P- = P+ and S^a = P^a + P^a + PL we finally get the simple 
result 

Nia = A+A'P+ + A-A'P- + A+A-Pl (S.19) 
The definition of the magnetic field can now be used to obtain the decomposition 

B,a = B+P+ + B-P- + B^Pl. (S.20) 



We get 



and 



B+ = i^A+ + '-A+ - '-A^ - A+^o, (5.21) 

dr r r 

B- = -i-^A- - -A- + - A-A"" (B.22) 
dr r r 



B° = -A+ - -A- - A+A-. (B.23) 



Now we include the trivial factor ^ in the definition of the vector potential using the 
notation A~^ = ^a~^, etc. We get 



B+ 
B- 

and 



l(i^a+--a\i + a+)], (5.24) 
r \ dr r ) 

- f-z4^a- + -a°(i-a-)^ (S.25) 
r \ dr r J 



B^ = ^{ia+ -ia- -a+a-). (5.26) 

One realizes at this point that the form of the ansatz (16) is simple in terms of the 
magnetic field giving 

a+ + i = ue''^, {B.27) 

a- -i = ue-'^, (S.28) 

a° = w. (5.29) 
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Using this we arrive at 



a"*" — a = 2iu sin </> — 2i, 



a'^a = + iue^^ — iue 



+ I = l- 2usva6 + u^ 



whence finally we obtain 



r \ dr dr r 



1 [ .du # 1 

B = - -I— e ^ - w— - e wue 

r \ dr dr r 



■i(f) 



and 



(5.30) 
(S.31) 

(S.32) 
(S.33) 

(EM) 



To obtain the magnetic fliix and the magnetic charge is easy observing that and P^^ 
are orthogonal to n^. So we arrive at eq. (22) 



(S.35) 



In order to calculate the energy density we note that since B = (-B+)*, P-^^ = (P^)*, 



\BiaBia = ■ B + • B^ 



is whence 



.du /d(f) 
dr \dr r 



u 



2\2 



Finally it gives rise to the following reduced dimensional action 



1 47r 



dr 



/ \ 2 / \ 2 

I du\ , 2 I d4> , w\ 1 



^dr ) i(ir~'"rj 2r^ ^ 



(S.36) 



(S.37) 



(S.38) 



as presented in eq. (24). 

It remains to calculate the action of an infinitesimal gauge transformation on the 
monopole ansatz. From the general expression 
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(5.39) 



we conclude that A5 = • A(r) is the only form not leaving the ansatz's configuration 
space, because only the derivative of nj is expressable through a combination of and 
P°-. Noting that 

Diab = Sabdi - eacbAic (5.40) 

and 

eacbrib = -iP^c + iPac (^-41) 

we obtain 

which, upon using the generic form (A.l) for the vectorpotential, yields 

5A+ = -A - iA+ ■ A, (S.43) 
r 

5 A- = -A + zyl- • A (BAA) 
r 

and 

5A^ = ^A. 
dr 

Replacing now the expressions for and A^ we finally arrive at 

6u = 0, (S.45) 

50 = -A (BAG) 

and 

6w = r-^A, (BA7) 
dr 

as claimed in eqs. (28). 
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Appendix C 



In this appendix we briefly describe the methods we employed to obtain the numerical 
results of this paper. 

The classical monopole-soliton as a solution of eq. (43) can not be obtained by direct 
numerical integration, because the fixed boundary condition u{oo) — 0,tt(0) = 1. Instead 
we applied the energy-functional relaxation updating the field tt„ = uo{x) known on a 
lattice of Xn = X ■ n ■ Ar = n ■ Ax points due to the conjugate gradient method 




(c.i) 



using e = 0.01, Ax = 0.01 - 0.05 and 200 - 1000 grid points. We obtained the same solution 
by a second order conjugate gradient method. 

We also found the eigenvalues of the effective potential (59) for monopole-like fiuctu- 
ations numerically. The first method employed a "shooting" algorithm, integrating from 
X = and from a; = oo up to a matching point, and modifying the initial derivatives in 
order to fit the logarithmic derivative at the matching. The second method directly located 
the zeros of the determinant of the discretized matrix 

^mn - <^'^Ax'^6mn = -Sn,m+1 " (^n,m-l + (^mn (2 + Ax"^ {Vefl[Un , u] - O;^)) (C.2) 

in the interval < a; < 1. The evaluation of the determinant made use of the simple 
recursion formula for tri-diagonal N x N matrices 

Dn = Qn ■ Dn-1 - Dn-2, (C-3) 
where qn is the A'"-dependent diagonal element and Dq = 1, Di = qi start the recurrence. 
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Figure Captions 

Fig. 1 The classical magnetic monopole soliton uo{r) which minimizes the energy func- 
tional E[u] + A£^[w] is shown as solid line. The resulting monopole charge con- 
tained inside a sphere of radius r is represented by the dashed line. 

Fig. 2 The effective potential for radial magnetic monopole fluctuations 5u{r) is plotted 
as a function of the scaled radius x — Xr. The ground state energy lOq — 0.95 in 
this potential is indicated by a horizontal straight line. 
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